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ABSTRACT 



Context. Optical modules for X-ray telescopes comprise several double-reflection mirrors operating in grazing incidence. The con- 
centration power of an optical module, which determines primarily the telescope's sensitivity, is in general expressed by its on-axis 
effective area as a function of the X-ray energy. Nevertheless, the effective area of X-ray mirrors in general decreases as the source 
moves off-axis, with a consequent loss of sensitivity. To make matters worse, the dense nesting of mirror shells in an optical module 
results in a mutual obstruction of their aperture when an astronomical source is off-axis, with a further effective area reduction. 
Aims. To ensure the performance of X-ray optics for new X-ray telescopes (like NuSTAR, NHXM, ASTRO-H, IXO), their design 
entails a detailed computation of the effective area over all the telescope's field of view. While the effective area of an X-ray mirror is 
easy to predict on-axis, the same task becomes more difficult for a source off-axis. It is therefore important to develop an appropriate 
formalism to reliably compute the off-axis effective area of a Wolter-I mirror, including the effect of obstructions. 
Methods. Most of collecting area simulation for X-ray optical modules has been so far performed along with numerical codes, involv- 
ing ray-tracing routines, very effective but in general complex, difficult to handle, time consuming and affected by statistical errors. 
In contrast, in a previous paper we approached this problem from an analytical viewpoint, to the end of simplifying and speeding up 
the prediction of the off-axis effective area of unobstructed X-ray mirrors with any reflective coating, including multilayers. 
Results. In this work we extend the analytical results obtained: we show that the analytical formula for the off-axis effective area 
can be inverted, and we expose in detail a novel analytical treatment of mutual shell obstruction in densely nested mirror assemblies, 
which reduces the ofl'-axis effective area computation to a simple integration. The results are in excellent agreement with the findings 
of a detailed ray-tracing routine. 
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1. Introduction 

Optical modules for X-ray telescopes consist of a number of 
grazing incidence mirror shells with a common axis and fo- 
cus. In a widespread design, the Wolter's, the mirrors comprise 
two consecutive segments, a paraboloid and a hyperboloid, in 
order to concentrate X-rays by means of a double reflection. 
The two reflections occur at the same incidence angle for X- 
rays coming from an on-axis source at astronomical distance 
(Van Speybroeck & Chase il9721) . The optical design of the mod- 
ule is primarily dependent on the required effective area, which 
determines the telescope's sensitivity. 

The most representative indicator of the concentration power 
of an optical module is assumed in general to be the effective 
area on-axis. However, the effective area in general decreases 
as the X-ray source moves off-axis. This in turn diminishes the 
telescope's sensitivity, hence may represent a severe limitation 
to its field of view. For this reason, the prediction of the on- and 
ofF-axis effective area modules is a very important task in the de- 
velopment of new X-ray telescopes such as NuSTAR (Hailey et 
al. |2010), NHXM (Basso et al. 2010), ASTRO-H (Kunieda et 
al. l2010j . and IXO (Bookbinder 2010). 

While the effective area on-axis is in general easy to calcu- 
late for a typical Wolter-I mirror module configuration, the com- 
putation becomes more difficult off-axis, because of the variable 
incidence angles over the two reflecting surfaces and the vari- 
able fraction of singly-reflected X-rays that are not focused and 



contribute to the stray light (see e.g., Cusumano et al. 120071 ). 
To make things worse, the assembled mirrors can shade each 
other if they are not spaced enough, which contributes to an even 
steeper decrease in the collecting area ofF-axis. A widespread 
method for performing the calculation, accounting for all these 
factors, has made use of accurate ray-tracing routines (see, e.g., 
Mangus & Underwood 119691 Zhao et al. 120041 ). which recon- 
struct the paths of a selection of X-rays impinging a mirror 
module. These numerical codes are in general accurate, be- 
cause they simulate the real incidence of rays on the optical sys- 
tem. However, they are complex and time consuming, especially 
whenever the simulation includes wide-band multilayer coatings 
to extend the reflectivity beyond 10 keV (Joensen et al. 119951 
Tawara et al. 1998 ). The reason is that the multilayer reflectivity 
computation, which has to be performed for every ray traced, is 
a complex procedure especially for wideband multilayers, which 
comprise many (~ 200) couples of layers. 

For this reason, although the ray-tracing approach should 
not be disregarded, it is interesting to derive analytical formu- 
lae for the effective area ofF-axis. In attempting to achieve this. 
Van Speybroeck & Chase (1972 ) discovered, by analyzing the 
results of ray-tracing simulations, that the geometric collecting 
area of a WoIter-I mirror decreases with the ofF-axis angle 6 of 
the source, according to the Formula 



AUO) = A»(0) 1 - — 
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where Aco(O) is the on-axis geometric area and ckq is the inci- 
dence angle for an astronomical source on-axis. 

More recently, a semi-analytical method for computing the 
off-axis effective area has been applied to solve the problem 
of optimizing multilayer recipes to the telescope's field of view 
(Mao et al. fT999l Mao et al. l2000l Madsen et al. l2009l l. In this 
approach, the effective area is computed by integrating the mir- 
ror reflectivity over the incidence angles, after weighting it over 
an appropriate function. Wine, derived from a ray- tracing. 

In a previous paper (Spiga et al. i2009l , we akeady derived 
a completely analytical method for computing the off-axis ef- 
fective area of a Wolter-I mirror with any reflective coating. In 
a subsequent article (Spiga & Cotroneo 2010), we developed 
this formalism, deriving the analytical expression of the afore- 
mentioned Wine function, which represents the distribution of 
the off-axis effective area over the incidence angles of rays, and 
we also used it to face the problem of multilayer optimization 
(Cotroneo et al. 1201 Oi l. The results were accurately verified as 
well, by comparison with the outcomes of a ray-tracing program. 
However, in these previous works, we did not consider the mu- 
tual shading (also known as vignetting) of mirrors, which may 
occur in mirror modules when shells are densely nested together. 
While in general the mirror module is designed to avoid any vi- 
gnetting on-axis, the problem may arise for sources off-axis and 
cause a further loss of effective area. Consequently, the results 
could be applied only to single mirror shells, or to mirror assem- 
blies for which the mutual obstruction is known to be negligible 
over all the field of view. 

In this work, we overcome these limitations and extend the 
developed formalism to the general case of obstructed Wolter-I 
mirrors in X-ray optical modules. We still assume that the mir- 
ror profile can be approximated by a double cone as far as the 
sole effective area is concerned, a condition in general fulfilled 
by optics with large /-numbers. In Sect. |2] we briefly review 
the results obtained for unobstructed single miiTor shells, and in 
addition we show how the analytical formalism can be inverted 
to derive the product of the two reflectivities from the desired 
effective area variation with the off-axis angle. In Sect. |3] we 
describe the geometrical parameters driving the nested mirror 
obstructions, derive the expression of the vignetting coeflicients, 
and obtain an analytical formula for the off-axis effective area of 
an obstructed mirror (Eq. (l38ll). We then derive in Sect.|4]some 
analytical expressions for the obstructed geometric area, and we 
apply the results to the geometrical optimization of the module. 
In Sect. |5] we prove the validity of the analytical formulae by 
means of a ray-tracing routine. The results are briefly summa- 
rized in Sect. |6] 

2. Analytical formulae for unobstructed Wolter-I 
mirrors 

We consider a Wolter-I mirror shell (Fig.lTJ with focal length /. 
We denote with Rq the radius of the reflective surface at the inter- 
section plane and with qto the incidence angle at the intersection 
plane for an astronomical source on-axis. They are related by the 
well-known relation 

/?o=/tan(4ao). (2) 

We hereafter assume that all incidence angles are shallow, there- 
fore tan(4Q'o) ~ 4ffo- The optical axis of the miiTor is oriented 
in the z-axis direction. We define L\ and L2 to be the lengths 
of the parabolic and hyperbolic segments (hereafter named pri- 
mary and secondary), which are supposed to be undeformed and 




/ 

to focal plane 



Fig. 1. Sketch of a grazing-incidence Wolter-I mirror shell with 
an off-axis X-ray source. We also show a ray direction vector 
before reflection (fe^), after the first reflection ik^), and after two 
reflections {k^. 

very smooth. We denote with Ry[ and the mirror radii at the 
entrance and exit pupil, respectively. We define o(a) to be the 
coating reflectivity for a generic incidence angle a, at the X-ray 
wavelength A, and assume that the geometric optics can be ap- 
plied, i.e., that A is small enough to avoid aperture diffraction 
effects but large enough for the scattering due to surface rough- 
ness to be negligible (Raimondi & Spiga 2010). The source is 
assumed to be at a finite, although very large, distance D » L\ 2- 
Finally, > Q denotes the off-axis angle, i.e., the angle formed 
by the source direction with the optical axis. The jc-axis of this 
reference frame is oriented such that the source lies in the xz 
plane, at jc > and z > 0. 

2.1. A review of previous results 

We briefly review the results about the computation of the off- 
axis effective area of unobstructed Wolter-I miiTor shells, derived 
in previous works (Spiga et al. i2009l Spiga & Cotroneo 1201 Oi l. 

1 . Owing to the shallow incidence angles, the double cone ap- 
proximation is in general applicable to compute the effective 
area of Wolter-I mirrors. The error in the effective area esti- 
mation, by replacing a Wolter-I profile with a double cone, 
is on the order of L/f or smaller, i.e., a few percent in real 
cases. 

2. For a source off-axis, the incidence angles on the primary 
and secondary segment, a[ and a2, essentially depend on the 
polar angles (p[ and (p2 of the points where the ray is reflected 
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(Fig. [T]). As long as is small, (pi ^ (p2- We thus denote 
their nearly-common value with (p. If the polar angle tp is 
measured from the ;c-axis, ai and 02 have the expressions 

ai = a()+6-6cos(p, (3) 

Q'2 = ao - 6 + 6 cos tp, (4) 

where 6 - Rq/D is the X-ray semi-divergence due to the 
distance of the source: for an astronomical source, 6-0. 
Equations © and dUl are valid if ai > 0, 02 > 0. We note 
that 



ai + 02 = 2ao. 



(5) 



3 . The geometric ratio of rays that undergo the double reflection 
to those impinging the primary segment is expressed by the 
vignetting factor 

Viw) - (6) 
Lia\ 

with the constraint that < Viip) < 1 . 

4. Using Eq. (|6]l, it is easy to derive an integral formula for the 
effective area of the mirror shell 



r {La. 
Jo 



)miarA(ai)ri(a2) dip. 



(7) 



where 

(Z^Qf) min — max{0, min[Liai(i^), L2<3'2(^P)]) ■ (8) 
5. In the frequent case Li = L2 (= L), Eq. (|7]i reduces to 



Ad{A, 0) = 



Jo 



I 



(9) 



where amin - min[ai(i^), Q'2(i^)] if positive, and zero other- 
wise. The expression of the o-min angle can be written in a 
compact form 

^min - max(0, ao -\6 - Ocos (p\), (10) 

In the ideal case of a perfectly-reflecting mirror, i.e., setting 
r = 1 for any a and A, the integral in Eq. (|9]l can be solved 
explicitly. For instance, if 6 - and = 0, we obtain 

AUO)^27iRoLao, (11) 

whereas, if 6 - and < 9 < ao, 



26 

A^(0) = Aoo(O) 1 1 

nao 



(12) 



where - and heretofore - we drop the dependence of the A's 
on A to denote the geometric mirror areas. Equation (fTSI) . 
after approximating n x 3, becomes Eq. ([T]), the expression 
found by Van Speybroeck & Chase in 1972 For 9 > ao, the 
geometric area has the more complicated expression 



A^(0) = A^Q) 



1-i 




n 


ao 



02 



ao 



(13) 



7. If (5 > 0, more analytical expressions can be obtained by 
solving the integral in Eq. The predictions are in very 
good agreement with the findings of detailed ray-tracing rou- 
tines. For more details, we refer to Spiga et al. 120091 

8. If 5 = and < < ao, we obtain an alternative form of 
Eq. (|9]), by changing the integration variable from 1^ to ai, 

A^(A,0) ^4RoL dai. (14) 

Jao-e _ (-Qjp _ Q,j-)2 

This equation can be used to derive the eff'ective area as a 
function of 0, at a fixed A. 
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Fig. 2. The product of the reflectivities of the two segments of a 
Wolter-I shell with ao = 0.212 deg, at 30 keV. The reflective 
surface is a multilayer coating, a Pt/C graded stack with 200 
couples of layers. The d-spacing in the multilayer follows the 
power-law model d(k) - a{b + ky^ (Joensen et al. 119951 ), with 
a = 77.4 A, b = -0.94, c = 0.223, a Pt thickness ratio r= 0.42, 
and a 4 A roughness rms. The solid line is directly computed 
from the multilayer structure, whilst the dots are the result of the 
inverse computation from the effective area in Fig. [3] 
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Fig. 3. The eff'ective area at 30 keV (A = 0.4 A) of a Wolter- 
I mirror shell with Ro - 148.5 mm, L - 300 mm, / = 10 m, 
ao - 0.212 deg, and a reflectivity product as shown in Fig. |2] 
(solid line). The effective area was used to re-derive the Tx(ai) 
function (Fig.|2] dots). 



2.2. Inverse computation: from Aoo(/l, 0) to the mirror 
reflectivity 

As a further development of the equations reported in the pre- 
vious section for unobstructed mirrors, we hereafter provide the 
inverse formula of Eq. (fT4l i. For a given A, from the effective area 
variation with in [0, ao], for 5 = 0, we derive the product of the 
reflectivities of the two segments 



def 



T^iai) = o(Q'i)r^(ff2), 



(15) 



where a2 - 2ao - ai (Eq. Owing to the symmetry of a\ 
and Q'2 with respect to the y-axis when i5 = 0, it is sufficient to 
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compute T^(ai) for < ai < cq. To this end, we define Ej^(ff) to 
be the ratio of the effective area to the on-axis geometric area 



Ea(0) = 



AUO) ' 



(16) 



and then it is easy to demonstrate (see Appendix |A| that the 
2"i(a^i) function, for ai G (0, ao), can be computed from the in- 
tegral equation 



ri(ai) 



— ( sin f 

ai Jo 



—(e-EA0)) 
a6 



dt, 

9=6(1) 



(17) 



where the expression in the [ ] brackets, after computing the 
derivative, is to be evaluated at 



6(t) - (ao - Q-Osinr, 



(18) 



where < f < n/2 is a dummy integration variable. The con- 
dition < cKi < ao then implies that < 6 < ao'. in other 
words, the £^(6) function has to be known in the entire interval 
(0, ao)- By computing the integral in Eq. ( fTTT i for any normalized 
effective area function E,i{6), one obtains the corresponding re- 
flectivity product at A (Eq. (fTsTi). for ai taking on values in the 
same interval. 

As a first example, we put in adimensional form the geomet- 
ric vignetting for a source at infinity, Eq. (fTZt . 



26 

E(0) = 1 , 

nao 

and substitute this into Eq. ( fTTl l. We obtain 



Txiai) 



ai Jo 



sinf - 



4( q-o - ax) 
nao 



sin^ f df, 



which can be immediately solved, yielding 
r,(ai) = 1. 



(19) 



(20) 



(21) 



This means that r^(Q') = 1 for any incidence angle and wave- 
length, as expected. 

As a second example, we show the computation of a reflec- 
tivity product from the effective area of a mirror shell with a mul- 
tilayer coating. The product of the two reflectivities at 30 keV, 
directly computed from the multilayer recipe, is shown in Fig.|2l 
as a function of a\. The effective area of the mirror shell at the 
same energy (Fig. [3J was computed in (0, ao), using Eq. (fT4b . 
Finally, we used Eq. (fTTb to re-derive the T,\{a\) function from 
the effective area curve. The result of the inverse computation 
(Fig. |2l dots) closely matches the original reflectance product 
(Fig.H line). 

3. Obstructions in nested mirror modules 

We now deal with a quantification of the obstructions that re- 
duce the off-axis effective area of a miiTor shell, when nested 
in a mirror module. Even if this effect is essentially a geometric 
shadowing, the effective area reduction depends on A, because 
the incidence angles, and consequently the mirror reflectivity 
ri(a), vary over the miiTor surface when the source is off-axis. 
This is especially true when multilayer coatings, which exhibit a 
complicated r^ia) function, are adopted. A reduction of the ef- 
fective area is then relevant for the X-ray wavelengths that were 
reflected at the obstructed regions. 



The clear aperture of a mirror shell can be obstructed by sev- 
eral factors: the dense packing of shells, the structures for me- 
chanical support of mirrors, and the presence of pre-collimators 
designed to reduce the stray light. However, we hereafter limit 
the discussion to the first point, i.e., the mutual obstruction of 
nested mirror shells, assuming that they are all co-axial, co-focal, 
and all have the same intersection plane at z = (Fig.|4]l. 

3. 1 . Obstruction parameters 

We now draw our attention to a particular mirror shell in the 
mirror module (Fig. |4]i, and adopt the same notation presented 
in Sect. |2] The obstruction of this shell (which we refer to as 
a reflective shell) can be assumed to be caused solely by the 
shadow cast by the adjacent shell with smaller radius, which 
we refer to as blocking shell. The radius of its outer (i.e., non- 
reflective) surface at z = is denoted with 7?^, which is forcedly 
smaller than Ro- We then define 7?^, R*^, L*, and L*^ to be respec- 
tively the maximum radius, the minimum radius, the primary 
segment length, and the secondary length of the outer surface 
of the blocking shell. We admit that, in general, L\ + L* and 
L2 + E*^. If Lj - E*^, we denote their common value with E* . 

A simple geometrical construction shows (Fig.|4| that an ob- 
struction of the reflective shell can occur in three ways: 

1 . If rays intersect the blocking shell's primary segment before 
impinging the reflective shell; 

2. If, after the first reflection, rays impact onto the outer surface 
of the blocking shell's primary segment; 

3. If, after the second reflection, rays are blocked by the outer 
surface of the blocking shell's secondary segment. 

The obstruction amount depends on the closeness of the two 
mirrors. This is often expressed along with FF, the filling factor 



FF = — . 



^0 



(22) 



The configuration with FF - 1 causes each shell to exactly fit 
the clear area section of the adjacent shell, hence it maximizes 
the effective area on-axis at the expense of the off-axis area. If 
FF > 1, the mirror assembly is self-obstructed even on-axis. For 
this reason and to enlarge the field of view of the optics, solutions 
with FF < 1 are in general envisaged. 

Nevertheless, it is convenient to adopt other parameters in 
the estimation of obstructions. The first of these is 



<D = 



Ro 



El 



+ ao, 



(23) 



which represents the maximum angle visible from the reflective 
shell at the intersection plane (see Fig. |5]l, through the entrance 
pupil. Clearly, <1> > o-o whenever < 1, and O < ao if FF > 1. 
As we later see, the <i> parameter drives the obstruction of the first 
kind, i.e., at the entrance pupil. 

The obstruction of the second kind, i.e., at the intersection 
plane, is chiefly determined by another parameter. 



Ro — Rn 



(24) 



which denotes the clear angular aperture at the intersection 
plane, as seen from the maximum diameter of the reflective mir- 
ror shell (see Fig. |5]l. The third obstruction parameter, which 
drives the third species of obstruction, is 



Ro - R*j, 

^ = — Tz 3ao, 

^2 



(25) 
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Fig. 4. Obstruction in an assembly of mirror shells. Rays im- 
pinging a mirror shell can be blocked by the adjacent shell with 
smaller diameter in only three ways, as listed in the text. The 
impact points are highlighted. Other mirror shells are not shown. 
The obscured regions of the primary mirror are indicated in gray: 
angles and the mirror spacings are greatly exaggerated. 

which represents the angular aperture visible from the reflective 
shell at the intersection plane (see Fig.|5]), through the exit pupil. 
The importance of these angles becomes clearer in Sect. 13.31 
when we derive the general formula for the obstructed mirror 
effective area (Eq. (l38ll). which depends on O, *F, and E as pa- 
rameters. 

The obstruction parameters can be related to each other 
O = + (ao - aa), (26) 



■ 3(ao - ^o), 



(27) 



where is the on-axis incidence angle on the blocking shell. 
Using Eq. (|2]l, these relations can also be written as 



(D 



-it- 



1 

4/ 

3 

4/ 



(28) 



(29) 





Fig. 5. Geometrical meaning of the obstruction parameters, O, 
*P, and 2, for a pair of Wolter-I nested shells. 



- if l: < L, 



<Ll, 



then S < 'P < O; 



- if = Li = LI and / » Li, then S >P ; 



3.2. Vignetting coefficients 

We now provide analytical expressions for the obstructions 
caused by the mirror nesting. To this end, we introduce a vi- 
gnetting coefficient, V„((p), with n = 1,2,3, for every kind of 
obstruction as listed in Sect. 13. II For an infinitesimal mirror sec- 
tor located at the polar angle ^, the «''' vignetting coefficient is 
defined as the fraction of primary segment's length left clear by 
the n''' obstruction. This definition is analogous to that of V{ip), 
the self-vignetting factor for double reflection (Eq. (|6]l), which 
we already treated in detail (Spiga et al. |2009l l. We then have, in 
addition to V{(p), three vignetting coeflicients: 

- Vi {(fi), due to the shadow cast by the blocking shell's primary 
segment before the first reflection; 

- V2{(p), due to the shadow cast by the blocking shell's primary 
segment after the first reflection; 

- V^iifi), due to the shadow cast by the blocking shell's sec- 
ondary segment after the second reflection. 

Firstly, we consider all vignetting factors to be independent of 
each other: the respective obstructed (i.e., lost) fractions of pri- 
mary mirror length are denoted with Q„ = 1 - V„. 

As we anticipated in the previous section, it is the O, *P, 
and E parameters that primarily determine the vignetting coefli- 
cients. More precisely, they are functions of the incidence angles 
ai and q'2, and they can be computed via the following formulae 



From Eqs. (1281 ) and ( l29b it follows, in particular, that 



Vi(ai) 
V2(ai) 



= 1 + 

= 1 -1- 



^;(o-gi) 

Liai 



L*(l, - a2) 



(30) 
(31) 
(32) 



6 



D. Spiga: Optics for X-ray telescopes: analytical treatment of the ofF-axis effective area 



> 0.4 











/ 


















/ 










> ^ ^ 










/ 




















/ 

/ 




















/ 




















/ 






































/ 






ubie reflection 
itrance pupil 
tersect. plane 








/ 






- V, dc 
■ V,, e 
■■ Vj, in 












/ 
















; 






Vg , UAH pupil 


























1 1 1 1 


, , , 




1 1 1 


1 1 1 


1 1 1 



20 40 60 80 100 120 140 160 
Polar angie (deg) 



Fig. 6. Vignetting coefficients as a function of the polar angle 
for L2 = 310 mm, L\ - 300 mm, - 290 mm, L\ - 280 mm, 
ao = 0.3 deg, 6^0 deg, Q = 0.25 deg, <D = 0.43 deg, = 0.4 deg, 
and S = 0.408 deg. We also plot V, the vignetting factor for dou- 
ble reflection (Eq. (|6]l). 
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Fig. 7. Vignetting coefficients as a function of the polar angle for 
the same configuration as in Fig.|6] but this time the source is at 
a finite distance {D - 127 m, 6 - 0.095 deg). We note that the 
obstruction at the exit pupil is less severe than for a source at 
infinity. 



where a\ +02 - 2ao, and with the usual constraint < V,, < 1 
for each «, otherwise we set V„ to or 1, respectively. The deriva- 
tion of Eqs. (l30t to (l32b is not difficult, but quite lengthy, so it 
has been postponed to AppendixiBl 

The dependence of the coefficients on (f is obtained by sub- 
stituting the expressions of ai and 02 (Eqs. (O and (|4|). We note 
that the obstructed region of the primary segment in the first and 
third kind of vignetting is located near the intersection plane, 
whereas V2 results from an obscuration of the primary mirror 
near z - +Li (see Fig.|4]i, in a precisely similar way to V. 

Example plots of V], V2, V3, and V, as functions of ip, are 
drawn in Figs. |6] and |2l The adopted values correspond to the 
case of two mirror shells with /?() - 210 mm, 7?^ - 207.9 mm, 
and / = 10 m: the other parameter values are reported in the 
figure caption, including the lengths of the mirrors, which have 
been chosen to fulfill the relations L\ < < Li < La- As we 
later see (Sect. 14.21 ). this choice is not accidental: it represents a 
compromise to minimize all obstructions. 

Equations ( l30t to ( |32] | and Figs.|6] and|2]show that: 

- If the conditions ai < <I), ffi < *f , and aa < ^ are fulfilled at 
a given ip, there is no obstruction at that polar angle; 

- If ai > *P, but y < Va as in Figs. |6] and |7] the obstruction 
at the intersection plane is ineffective, because all rays that 
were blocked would have missed the secondary segment; 

- The obstructions related to V\ and V2 are maximum at polar 
angles close to n; 

- The effect of Vt, is larger at ^ » 0, where a\ is shallower and 
a2 is larger; 

- With the source at infinity, the obstruction related to can 
be very large, especially near (/? a; 0. In contrast, the effect of 
V'3 is mitigated if the source is at a finite distance, because 
a2 becomes smaller; 

- If Li > Lj, then V\ > V2 for all ip; 

-If Lx = L\ and / » Li, then ^ k (Eq. ^) and also 
V\ a V2 for all If. 



3.3. General formula for the effective area of obstructed 
mirror shells 

We already pointed out that the obstructions described by V and 
V2 are located close to the maximum diameter of a Wolter-I mir- 
ror shell, whereas the ones related to Vi and V3 concern mainly 



the region close to its intersection plane. Therefore, V and V2 
are in competition for the obstruction near z = +Li , while Vi 
and Vi do the same near z = (Fig.|4ll. The total obstruction at 
the generic polar angle if is then 



gtot = max(e, Q2, 0) + max(Qu Q3, 0), 



(33) 



where the "0" value has been added to avoid negative obstruc- 
tions. The corresponding total vignetting factor is 1 - Qtot, i-C-, 



Vtot = min(y, V2, 1) - max(l - Vj, 1 - ^3,0). 



(34) 



Using the expressions of V (Eq. (|6]l) and V2 (Eq. (ISTTi). we can 
write the first term of Eq. ( |34] | as 



(1202 *!' 
min(y, V2, 1) = mini , — , 1 



Using Eqs. (l30b and (l32l l. the second term reads 



(35) 



(LUai -O) L:(q'2 -E) \ 
max(l -yi,l -y3,0) = max|^-;^ , ,0 .(36) 



Li a. 



If positive and not larger than 1, the total vignetting factor 
(Eq. (l34li). can be used to compute the obstructed mirror effec- 
tive area. An infinitesimal sector of the primary mirror with po- 
lar aperture A^, as seen from the off-axis source, has a vignetted 
geometric area RoVtotLiaiAtp, hence the total effective area is 



Ad(A, 0) = 2Rt 



f ViotLiai 
Jo 



rA(ai)rA(a2)d(p. 



(37) 



By substituting Eq. ( l34l l into Eq. ( l37b . it is straightforward to 
derive the final expression for the off-axis effective area of the 
obstructed mirror shell 

AoiA, 0) = 2Ro f [(La)™„ - e„,J>o r^ai) ^(az) d^, (38) 
Jo 

where 

(La)mm = min(Liai,L2Q'2,Li*I'), (39) 
gn^ax = max [L\(ai - (D), L;(a2 - S), 0] , (40) 

provided that > 0, a2 > 0, as usual. The []>o brackets in 
Eq. ( l38T l mean that the enclosed expression (La)niin - Qmax, if 
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Fig. 8. There is no vignetting at the intersection plane if the dis- 
tance of the ray reflected at the maximum and the minimum di- 
ameter is smaller than the spacing of the two mirror shells. 



This result has an immediate geometric interpretation, by noting 
that left-hand of Eq. (l44l l is exactly the maximum possible dis- 
tance from the reflective shell of a ray reflected twice (Fig. [8]l. 
All other rays undergoing a double reflection cannot exceed this 
distance, which therefore represents the minimum separation for 
the two shells at z = 0. 

Finally, we define L to be the equivalent length of the mirror 

2 1 1 

^ = — + — , (45) 
L Li L2 

and using Eqs. Q and (|4|, we can express the obstruction avoid- 
ance conditions as 

ao + 6 + < (46) 
ao-6 + 6<I., (47) 
aoL<Ro-Rl. (48) 

We note that if L\ - L2 - L, then also L - L, and that the last 
condition solely depends on the mirror pair geometry, not on the 
off-axis angle. 



negative, must be set to zero. For this reason, the integration can- 
not be carried out independently for the two terms. Moreover, 
because of the presence of ^' in Eq. (|39] ) and the different values 
of O and E in Eq. ( l40b . the computation is asymmetric with re- 
spect to the y-axis, also when 6-0. Therefore, the integral in 
Eq. dSST l is not equivalent to twice the same integral over [0, 7t/2], 
unless all these conditions are fulfilled (as in Sect. 14. II ): i) 6 - 0, 
ii) <£ ^ 2, and iii) there is no vignetting of the second kind. 

3.4. Conditions for obstruction-free mirror she//s 

We now discuss the conditions for the mirror shell to be 
obstruction-free. It can be seen, from Eqs. (39[ and (l40t . that 
there is no obstruction if, and only if, for all if 



min(Liafi, Laca) < ii*!*. 



(41) 



If these inequalities are fulfilled, Eq. (I38l l correctly reduces to 
Eq. O. The first two conditions are simply met if the maximum 
values of ai and 02 do not exceed O and Z, respectively. The 
third condition requires that either ai < ^ or a2 < jj-^, which 
in turn are equivalent to 



ai + 02 < 



i.e., using Eq. (|5]i, 

2aoL2 



Li +L2 



(42) 



(43) 



which reduces to ao < if Li - L2. The definition of *P then 
allows us to write Eq. ( l43l l as 



2L1L2 
Li +L2 



< Rn — Ra 



(44) 



' From Eq. d42b . it might alternatively follow that a\ < or 
Qr2 < ^- However, if Li < Li this would still imply that ai < 
whereas if L\ > Li we would still have 0-2 < T < In both cases, 
Eq. ( 142b remains valid. 



4. Some applications 

4.1. Analytical expressions for the geometric area 

As a first application, we derive some expressions for the ge- 
ometric area (i.e., in the ideal case ri(a) - 1 for all a) of an 
obstructed mirror as a function of 0. For simplicity, we only con- 
sider the case that <5 = and Li - L2 - L\ - L*2, and we 
suppose the mirror to be unobstructed on-axis, i.e., that <S> > oq 
(Eq. (I23I1). Finally, we reasonably assume that / » L, so that 
O X! X! Z (Sect. 13. Il l: consequently, > ao and there is no 
vignetting at the intersection plane (as in Fig.[8]l. Equation ( [39] ) 
then becomes, as in the unobstructed case. 



(La)n,in = L ■ min(ai, ai)- 



(49) 



Adopting O as a unique obstruction parameter, and since 
min(Q'i, a2) + max(Q'i, a2) - 2q'o (Eq. (|5])), Eq. ( l40t turns into 



gmax - L [2aQ - min(ai, ai) - <!>] : 



(50) 



on the condition that both terms are non-negative. Substituting 
this expression into Eq. ( l38l l and using Eq. (fTOl l with 5 = 0, we 
obtain 



Ad(0) 



= 2RoL [ [. 
Jo 



ffmin -/3max]>0 dip. 



where we defined 

Q'min = [ao -6\C0S(p\]>Q, 
jSmax = [ao+e\ COS ip\ - <^]>o, 



(51) 



(52) 
(53) 



and where the []>o brackets have the same meaning as those 
appearing in Eq. ( l38l l. 

We firstly consider the case ao < O < 2q'o, which implies 
that O - ao < ^o- As long as < <I) - ao, amin > but ySmax = 
for all (fi, i.e., the mirror is not obstructed and Eq. (fT2l) remains 
valid. 

We now increase 6. Since (D - ao < 0/2 < ao by hypothesis, 
we are allowed to consider the case O - ao < 6 < <I>/2. For 
all ip, we still have amin > 0, but this time /Smux is positive for 
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Fig. 9. Normalized geometric area, Aoo(0)/Aco(O) of an ob- 
structed mirror shell with ao - 0.2 deg, 6 = 0, and Li = La, as a 
function of the off-axis angle, for different values of the obstruc- 
tion parameter O, with O x ^' x 1, and L = L*. The curves are 
traced using the analytical formulae reported in Sect. 14.11 The 
curve for the unobstructed mirror (the solid line) is also valid for 
any obstructed mirror with (5 > 2ao. 



1 cos (p\ > <I> - ao: for these values of if, the mirror starts to be 
obstructed. Then Eq. (BTI ) changes into 



Aoc(0) = 4-R()L I (a^) - cos w) dw 

«yarccos — ^ 



pa] 

4RoL 

Jo 



(C) - 20 cos ip) dtp. 



(54) 



where the two terms are the unobstructed and the obstructed part 
of the area, respectively. We note that the two integrands are pos- 
itive in the respective integration ranges, hence the integration 
can be carried out and we obtain 



A„(0) 



= 1 - 



20 



1 + S 



(55) 



where we normalized the result to Aoo(O), the on-axis geometric 
area (Eq. (fTTT i). and defined the non-negative S (x) function. 



S (x) = Vl - x^ - X arccos x, with < x < 1 . 



(56) 



We recognize in the two first terms of Eq. ( fSSl l the usual unob- 
structed geometric area trend (Eq. (fT2li). whereas the term pro- 
portional to the S function represents the area lost because of 
the obstruction. Since 5(1) = 0, Eq. (l55l l converges to the unob- 
structed area trend at = O - q-q- 

We now consider the case > 0/2, still for <1> < 2q'() and 
(5 = 0. The first term of Eq. (|54] | remains unchanged, even if the 
integration range shrinks because of the larger values of 0. In the 
second integral, indeed, we have to replace the lower integration 
limit with arccos^ to keep the integrand non-negative. This in- 
tegration limit also guarantees that o-min > even for > ao- 
Because of this change, the resulting expression 



Aue) ^ ^ 20 

Aoo(O) nao 



l^S\^]-2sl^^ 



(57) 



has an additional term with respect to Eq. ( fSSl l. 



Finally, we assume that O > 2ao- This time < <1> - ao, 
hence the geometric area decrease deviates from linearity (for 
> ao) before the mirror begins to be obstructed (0 > O - ao). 
The condition ff^in > then implies that the lower integra- 
tion limit in the first term of Eq. ( |54| ) has to be replaced with 
arccos '-j-: consequently, a,^in = for < ^ < arccos ^, and 
the obstructed term is zero. We conclude that for <t) > 2ao and 
6 = the geometric area trend equals the unobstructed one, i.e., 
Eq. ( fT2l i for < ao, and Eq. ( fT3l l for > ao- As expected, in the 
limit O = 2ao Eq. ( l57b becomes identical to Eq. ( fT3b . 

As an example, we trace in Fig. |9] some geometric area 
curves of a mirror shell for different values of the obstruction pa- 
rameter O > ao, using Eqs. (fT2l i. ( l55b . and i5% in the respective 
intervals of validity. We also plotted the unobstructed geometric 
area of the mirror (Eqs. ( fT2l i and (fT3])), which is also valid for 
O > 2ao. Some curves are also validated in Sect.|5]by means of 
an accurate ray-tracing computation. 



4.2. Design of the mirror module 

The obtained results can also be applied to the problem of de- 
signing a mirror module. In general, the radius and the length of 
the outermost shell of the module is assigned on the basis of the 
allocable space for the optics payload. Starting from this one, 
shells with decreasing radii are added, leaving a sufficient spac- 
ing to minimize the mirror vignetting for off-axis angles within 
the field of view of the optic. It is then convenient to reduce the 
obstruction, by increasing O and 2 (Sect. 13.1b for every pair of 
consecutive shells. For any choice of Ro, Rq, and Li, this can be 
obtained by designing each mirror pair with L* < L\ (Eq. (|26]) ) 
and < Li (Eq. ( l27b ). On the other side, mitigation of the mir- 
ror self-vignetting for double reflection (Eq. (|6|) requires that 
L2 > Li and L*^ > L*. Hence, if we label the mirror shells with 
k = 1,2,... from large to small diameters, a performing module 
design might consist of mirror shells with 



< L 



(58) 



i.e., decreasing lengths as the diameter is reduced, and with sec- 
ondary segments longer than the respective primary but shorter 
than the primary of the adjacent mirror shell with larger diame- 
ter A module design with mirror lengths scaled to the diameter 
has already been studied by Conconi et al. (120101) to minimize 
the defocusing due to the field curvature in the WFXT telescope. 

Therefore, after choosing a sequence of mirror lengths ac- 
cording Eq. ( 1581 ). we apply the obstruction-free conditions 
(Sect.[I!4b. Using Eqs. and (|29]l, we rewrite Eq. (|47]i as 

ao,A +0< (Ro,k - Rom 1 - 1 ) ( - | , (59) 

where Tk+\ is the thickness of the (fe-n 1)''' shell, in general chosen 
to be proportional to the radius to maintain a constant mirror 
stiffness throughout the module. If we now define l/L2/,/(:+i = 
1 IL'2,k+\ ~ ^l^fi we can write the last equation as 

RoM\ +Tk+i < Ro,k - {ao,k + 0)L2f\k+i- (60) 
In reality, Eq. dSST l implies that Z < C) for all pairs of mirrors, 
hence Eq. ( l47T i also fulfills Eq. ( |46] |. The last relation to be satis- 
fied is Eq. ( 1481 ). which simply reads 

Ro,k+l +Tk+l ^ Ro,k - O!0,kLk+\. (61) 

We then derive from Eqs. (I6OI 1 and (1611 1 the condition to be ful- 
filled by the A;''' couple of shells, in order to avoid obstructions 
up to an off-axis angle 



RoMi + Tk+i < Ro,k - max \ao,kLk+uiao,k + 6')i2/,t+i] ■ 



(62) 
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Fig. 10. Initial positions and destinations of 40000 rays at the 
entrance pupil (points) for a Wolter-I mirror shell with L - L*, 
and the angular parameters reported in the legend (cD *P S). 
The radial scale has been expanded. Only rays that would have 
struck the primary mirror were traced. The limits of the regions 
of different vignetting (dashes) are computed from the vignetting 
coefficients. 
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Fig. 11. Comparison between some analytical curves of Fig. |9] 
(lines), and the results of an accurate ray-tracing (symbols) for 
an obstructed Wolter-I mirror with the same ao = 0.2 deg, 
Ro - 139.6 mm, L = L* - 300 mm. The outer radius of block- 
ing shell, takes on the values 137.24, 138.29, and 138.55 mm, 
for O to match the values used to draw the analytical curves (but 
the unobstructed one). The accord between the formulae and the 
tracing is excellent within the statistical error. 



If mirror lengths are chosen according to Eq. ( l58T l. the last re- 
lation enables the computation of the maximum possible value 
of Rq^u+x from /?o,i:, i-i.t+i, L2,k+i, and the relation = T{R()^k)- 
When applied recursively from the outermost radius inwards, it 
provides us with the optimal population of mirror shells in the 
optical module. 

5. Validation with ray-tracing results 

We hereafter validate the formulae found in the previous sections 
by comparing them with the results of a detailed ray-tracing. As 
a first example, we validate the expressions of the vignetting 
coefficients (Sect. I3.2l i. In Fig. [TO] we show the entrance sec- 
tion of an obstructed Wolter-I mirror shell with Rq - 139.6 mm, 
ao - 0.2 deg, / = 10 m, and L[ - L2 - 300 mm. Rays 
coming from a source at infinite distance, off-axis by an angle 
9 - 0.15 deg, have been traced by simulating their reflection 
on the mirror: the positions at the entrance pupil of rays that 
impinged the primary segment are drawn, disregarding the oth- 
ers. The points then fill the geometric section of the primary, as 
viewed from the direction of the source. The blocking shell has 
the same focal and length of the reflective one, but a different 
radius R^ - 138.3 mm. 

Depending on their initial coordinates, the traced rays can 
obstructed in various ways, highlighted with different colors in 
Fig. (To] or reach the focal plane if they fall in the red region of 
the mirror aperture. We note that, in agreement with the discus- 
sion in Sect. 13.21 the obstructions at the entrance (yellow) and 
the exit pupil (blue) are mainly located close to the intersection 
plane, i.e., the inner contour of the colored area, while the ob- 
struction after the first reflection (orange), and the vignetting for 
double reflection (green), mainly occur for rays firstly reflected 
at locations far from the intersection plane. Moreover, in this 
case (Y ^ O > ffo), the obstruction after the first reflection is 
ineffective because the orange region, which encloses the rays 
blocked after the first reflection, is completely surrounded by 



the green region, which comprises the rays that missed the sec- 
ond reflection. This is in more than qualitative agreement with 
the analytical tractation: in the same figure we have also traced 
the analytical expressions of the vignetting coefficients (dashed 
lines), after translating them into expressions of the radial co- 
ordinate (see appendix IbJ and projecting them at the mirror's 
entrance plane in the direction of the incident rays. We note that 
the lines perfectly follow the boundaries of the vignetted regions: 
the agreement shows that the expressions of the vignetting coef- 
ficients (Eqs. ( l30l l through ( l32l i) are correct. 

As a second example, we compare the analytical expressions 
of the geometric area found in Sect. 14.11 with the findings of 
the ray- tracing. Figure [TT] reports some of the normalized geo- 
metric area curves (lines) akeady shown in Fig. |9l i.e., for no 
obstruction and for two values of O in the interval [ao,2Q'o). 
We have also plotted as symbols the geometric area values ob- 
tained by applying a ray-tracing routine for a Wolter-I, reflective 
mirror shell and a co-focal blocking shell with the same length 
but variable values of /?y < Rq. Two values of R*q (138.55 and 
138.29 mm) were chosen to match the two values of O (0.2 and 
0.25 deg) we used to trace the analytical curves, while the small- 
est one (137.24 mm) was selected to return - 0.45 deg, i.e., a 
value larger than 2q'o- 

Inspection of Fig. [TT] shows that the findings of the two 
methods are in excellent agreement, within the error bars of the 
ray-tracing. On the other hand, the ray-tracing routine used is 
a quite complex code, and the relative computation required a 
few hours time to reach a few percent accuracy, whilst the ana- 
lytical curves can be traced immediately and without being af- 
fected by statistical errors. We also note that the results of the 
ray-tracing for O = 0.45 deg are perfectly reproduced by the 
unobstructed mirror analytical curve, in agreement with the con- 
clusion in Sect. 14.11 that there is no obstruction at any off-axis 
angle if O > 2Qfo. This happens because, with such a loose mir- 
ror nesting, 6 must be larger than oq for the mirror to start be- 
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Fig. 12. Effective area of an obstructed mirror shell with Rq - 
148.5 mm, / = 10 m, ao = 0.212 deg, L =L* =300 mm, and the 
X-ray source at infinite distance, 6 arcmin off-axis. The curves 
are computed for three different values of R^. The reflective shell 
has a Pt/C graded multilayer, whose d-spacing variation in depth 
is a power law, d(k) - a{b + ky, where k - 1 ... 200, a - 
1 15.5 A, = 0.9, c - 0.27. We moreover assumed a Pt thickness- 
to-d-spacing ratio F = 0.35, a roughness rms cr = 4 A, and a stack 
that ends with Pt on top. 



ing obstructed. In these conditions, the obstruction is ineffective 
because it is completely superseded at all polar angles by the 
vignetting for double reflection. 

Finally, we directly compare in Fig. [12] some obstructed ef- 
fective area off-axis curves, as a function of the X-ray energy, as 
computed from a ray-tracing and using Eq. (l38T l. The reflective 
mirror shell has a fixed radius and a graded multilayer coating 
(see e.g., Joensen 1995) to enhance its reflectivity in hard X- 
rays up to 50 keV and beyond. The characteristics of the reflec- 
tive shell and the description of the multilayer stack are reported 
in the figure caption. The off'-axis angle 9-6 arcmin is the same 
for all curves, while the outer radius of the blocking shell, R^, 
has been varied. As expected, the effective area decreases as the 
tightness of the nesting increases (i.e., as Rq approaches Ro), and 
moreover the findings of the ray-tracing (symbols) are in perfect 
agreement with those of the analytical computation (lines). 

We note that the high-energy part of the effective area is less 
affected by obstruction effects, because it results from reflection 
at polar angles (f x +7t/2, where max(Q'i , a2) takes on the small- 
est value, i.e., ao, and the reflectivity is higher: this is also the 
angular region where the obstruction is lower (see Fig.fTOli. 

6. Conclusions 

We have developed the analytical formalism for the off-axis 
effective area of Wolter-1 mirror shells, in double cone ap- 
proximation, which we began to describe in a previous paper 
(Spiga etal. ,2009j . 

We have shown that the analytical expression of the effective 
area off-axis can be inverted to derive the product of the reflec- 
tivity of the two segments (Sect. 12.21 ). This might be useful to 
future developments for computing a suitable multilayer recipe 
to return the desired effective area trend off-axis. 



We have found analytical expressions for the vignetting co- 
efficients (Sect. l3^ . for the three possible sources of obstruction 
in nested mirror modules, as a function of the azimuthal coordi- 
nate of the mirror surface. 

Using the vignetting coefficients, we have derived an integral 
formula (Eq. ( l38l l) for the obstructed effective area of a Wolter-I 
X-ray mirror in double cone approximation, with any reflective 
coating, including multilayers. The computation only requires 
the standard routines for the reflectivity of the coating, and an 
integration over the azimuthal coordinate of the mirror shell. 

We have obtained analytical expressions of the obstructed 
geometric area (Sect. [431 1 for the case of a source at infinite dis- 
tance, and applied the result to the problem of designing an op- 
tical module that does not suffer from the mutual obstructions of 
mirrors (Sect. [431 1. 

Finally, the results have been validated by means of a com- 
parison with the findings of a detailed ray-tracing (Sect.|5]). 

As a final application, we note that each vignetting coeffi- 
cient can be adapted to estimate the unwanted vignetting caused 
by collimators aimed at reducing the stray light in mirror mod- 
ules (see, e.g., Cusumano et al. 2007). For example, V] would 
quantify the vignetting of the baffle at the entrance pupil, if R^^ 
is interpreted as the outer radius of the collimator ring and L* as 
its distance from the intersection plane. With analogous substi- 
tutions, V3 would represent the vignetting of a baffle located at 
the exit pupil, and V2 would express the vignetting of the baffle 
at the intersection plane, even though this kind of baffle can be 
designed to avoid any obstruction of focused rays (Sect. [X4] l. 
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Appendix A: Inversion of the effective area integral 

We hereafter report the derivation of the inverse integral equa- 
tion (Eq. (fTTll). which allows us to derive the reflectivity product 
from the off-axis effective area (Sect. l2!2] l. We start from the ex- 
pression of the unobstructed effective area for Li - L2, 6 = 0, 
and < < ao (Eq. (fT4l i). which we rewrite in terms of the 
adimensional ratio E^(8) (Eq. (fT6]l). 

EA(e)^ , ^ ==dai, (A.l) 

where < < oq. By setting the positive variable co = oq - ai, 
the integral becomes 

2 r«(ao-w)7iM. 

Ea(0) = , doj. (A.l) 

^"Q Jo Ve^ - ^2 

We now set Ka{cij) = w^(q'o - (^)T,i(cj), u - l/a> and w - 1/0. 
Equation (IA.2I) can be rewritten as 

2.^0^=4 r^md. (A.3) 

The right-hand side of Eq. ( IA.3l l is the well-known Abel integral. 
It can be thereby solved for (see e.g., Stover[T995) 

Then, restoring the T function and the and co variables, we 
obtain 

Ta((o) = r f ^[0-EAmd0. (A.5) 

oj{ao - u) Jo Vw2 - 6l2 d0 
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Finally, we substitute w - co-ai in Eq. dA.Sb with < ffi < ao, 
and change the integration variable from to t, where = (ao - 
ai) sin /. Hence, we obtain the final result 



I sin t 

Jo 



d6 



(0 ■ EAO)) 



Moreover, we know from Sect. l2.1l also that (pi ^ (p2, thus we de- 
note with (f the nearly-common value of all these polar angles. 
We then rewrite Eq. ( IB .41 ). using Eq. ( IB. lb . as 

.(t - (ao + 6-ecosip) 



dt, 

0=0(1) 



(A.6) 



(B.6) 



where the [ ] brackets mean that the enclosed expression has to 
be evaluated at (ao - ai ) sin t. We have so obtained Eq. ( fTTb . 



Appendix B: Derivation of the vignetting 
coefficients 

We hereafter report the detailed derivation of the vignetting co- 
efficients (Sect. 13.21 ). The geometry of mirror shell obstructions 
and the meaning of symbols are explained in Figs.[T]|4] and|5] 

B. 1. Vignetting at tlie entrance pupii, Vi 

This obstruction is caused by the shadow cast by the primary 
segment of the blocking shell onto the primary segment of the 
reflective shell. For this to occur, the shadow of the edge of the 
blocking shell at - {R^cos(fo,R^sm(p(),L*^) has to intersect 
the reflective shell at some = (ri cos(fi,ri sini^i,Zi), with 
< Zi < Li . The shaded part of the mirror length is then Zi /Li . 

To find Z[((pi), we consider the generic ray emerging from 
the source, which is located atS_- [D sin 0, 0, D cos 6), and pass- 
ing by r^. Now, Tj, the ray intersection point on the reflecting 
primary, whose equation is 



ao + 6* - 6 cos if 

Clearly, Zi > only if ^ > 0. If the spacing of the two mirrors 
is small enough, we may also approximate 6* ~ 6. The first vi- 
gnetting factor is thereby Vi - 1 "I*-, i.e., recalling the definition 
ofai (Eq.dS), 

L*((D-ai) 



yi(ai) - 1 + 



L\ai 



(B.7) 



if positive and less than 1. We have so obtained Eq. ( l30l l. As 
usual, if this expression returns a negative value at some ip' , then 
V\(ip') - 0, or, if larger than one, V\(ip') - 1. If Li = Lf we find 
that V\ takes the simple form of 



Vi(ai) - — , 



(B.8) 



where in this case we also note that, if Vi were positive and 
smaller than one for all ip, the geometric area of the obstructed 
primary segment would become 



Ax(e)^ 2RoL 



ai 

Jo 



Vi(ai)dip, 



(B.9) 



Zi = , with ri > Ro, 

ao 



(B.l) 



in conical approximation, is identified via the vector equation 



(B.2) 



where x denotes a cross product. This returns two independent 
scalar equations 

(l - D ) '"1 '^1 + ^-^^^('"1 - ''^o) = Rm cos <Po - L\6 

\ . sin (On y ^ ^ . ' V ■ / 



(1-5) 



(n - Ro) =;?;jSin(^o 



where we defined 6* - R^^ /D to be the beam divergence at the 
blocking mirror shell. 

The impact position is obtained by solving the previous 
equations for ri and tfi. If = 0, the solution is expected to 
be independent of i^o, and the first points to be shaded are near 
the intersection plane, i.e., ri ^ Ro- To treat the general case of 
> 0, we search for a perturbative solution in the approxima- 
tion of small 6, so we set ri = Rq + ^ and (fi = (fo + s, with 
< ^ <*c /?() and |e| (fo- By substituting in Eqs. ( IB.3b . approx- 
imating to the first order, and solving, we obtain a linear system 
whose solutions are 



which - as expected - returns the area of the corona delimited 
by the largest radii of the two shells 

Am = InRoU^ = 2nRo{RM - R'm), (B.IO) 
where we have used the relation Rm - Ro + co^^i ■ 

B.2. Vignetting at fhe intersection piane, V2 

In this case, vignetting occurs after the first reflection from 
the primary segment of the inner shell, i.e., a point of the 
blocking shell's outer surface at z - 0, with coordinates 
£,) - (/?Q cos ^0,^0 sin 1^0,0), may intercept a ray reflected 
by the primary segment of the reflective shell at = 
(n cos^i,ri sin^i,Zi), with < Zi < Li. If it does, obstruc- 
tion occurs at Z\{(fi) < z < L\. The coordinate Zi is located via 
the vector equation 

x(£i -£()) = 0' (B-11) 
where is the ray direction after the first reflection. This unit 
vector has the following expressiorQ (Spiga et al. 12009b 

9 + {6 - 2ai)cos^i ' 

(B.12) 



(6- 



2ai) s\nip\ 
-1 



<^-{ao + 5-ecosipo) 
ao + 5* - 6 cos ipo 

£ — 

Ro (ao + 6* - 6 cos (po) 



'Sin (fo. 



(B.4) 
(B.5) 



where we have used the definition of <I> (Eq. ( l23T l) and neglected 
terms proportional to ao6. 

We note that s is of the order of or less, i.e., it is negligi- 
ble with respect to ipo itself. We can then assume that (pi » (fo. 



^ From the expression of the k^ vector, it is easy to derive the stray 
light pattern on a detector at a generic distance d, i.e., at z = -d. The 
final coordinates of rays generated close to the intersection plane are 

x{ip) = i?ocos(/j —d, 

= Ro sin (p---^d, 

and, substituting the expression of the vector components, we obtain 
the parametric equation of the pattern 

x{ip) = [Ro - (2ao + 6)d] cos ip + Qd cos lip, 
y(ip) = [Ro - (2ao + S)d] sin ip + dd sin 2ip. 
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After some manipulation of Eq. dB.llb . we obtain the two scalar 
equations 



(B.13) 



f (ffo cos ifx + kxx)r\ - kuRo + aoR*o cos ipo 
\ (ai) sin If I + kiy)ri = ku.Ro + aoRlsintpo ' 

and, using Eq. ( IB.12l i. we find that the solution of Eqs. ( IB.13I I 
for = is = (fii and 



rf' - R, + 



ao 



ao + 6 



(Ro-R'a). 



(B.14) 



Since Rq > we have rf' > Rq, hence zf' > (Eq. dRB ). 

If > 0, we set ri = rf^ + ^, ipo = ifi + s, and proceed as in 
Sect. IB. II The solution to a first order approximation is 

aniai) + 6 — at) 

^ " ^^'^0 - RD- (B. 15) 

We are not interested in the exact expression of s, which is of the 
order of 6, so we can assume again that ipo x ipi and neglect the 
i^'s subscript. In contrast, from ^ and Eq. ( IB. 14b we can derive 
an expression for ri , and using Eq. ( IB. II ) we obtain Zi 

^0 - ^0 



(B.16) 



which is always non-negative. All points with z > Zi are then 
obstructed. The resulting vignetting coefficient is V2 - Zi/Li, 
i.e., using the definition of (Eq. (l24l) ). 

Viiai) - — . (B.17) 
We have so obtained Eq. (ISTT i. 



B.3. Vignetting at tine exit pupii, Vj 

In this case, the obscuration occurs after the second reflec- 
tion, on the secondary segment of the blocking shell. A 
generic point of the blocking shell at z = -Lj, i e., - 
(R*^cos(po,R*^sinipo,-L2), may intercept a ray after it was re- 
flected at = ('"2 cos (p2, r2 sin (^2, Z2), with -L2 < Z2 < 0. If it 
does, the point £2 is located along with the usual vector equation 

K2 ^ d-i ' 



(B.18) 



where the z-coordinate of the secondary segment is given by 
'"2-^0 



Z2 = 



3ao 



-, with r2 < Rq, 



(B.19) 



and ^2 is the direction of the ray after the second reflection, 

-2 ~ -1 ~ ■ IkdlLi-' (B.20) 
because the tangential component is conserved, while the normal 
component to the surface reverses its sign in the reflection. To the 
small angle approximation, since (f2 - (fi (Sect. lZTT ). 



- cos ip\ ^ 

- smifx 
3ao 



(B.21) 



and fcj is given by Eq. ( IB. 12b . The final ray directiorQis then 

(B.22) 



k,^ 



-0 + (6 - 4ao) cos if I 
(5 - 4-ao) sin(/)i 
-1 



^ From Eq. l lB.22b we obtain that all rays reflected twice at loca- 
tions close to the intersection plane (r = Rq) converge to a single point 
F = (-Of, 0, -f), as expected, with 

^0 



/' 



4ao -(5' 



Setting (f2 - (fx in Eq. (IB. 18b . we obtain the two independent 
equations 



(Sao cos (fx + k2x)r2 = ^2.1-^™ + 3ao^m cos (fo 
(3ax) sin (fx + fe2v>2 = hyR'm + So-o^m sin (fo ' 



(B.23) 



where we set R[^ ^ /?o - 3ci'()L2- The solution for 0-Qis(fo = ifx 
and 



3ao 



(B.24) 



ao-o 

If > 0, we set r2 - r^"' + ^, (fo - (fx+sas in Sects. IB.ll 
and IB. 21 Substituting this into Eq. ( IB. 231 ), approximating at the 
first order and solving the linear system, we obtain 1^0 ~ (f\ and 
the expression of f . Adding this to rf' we obtain 



r2^R'^ + —(R'^-R*J, 

ai 



(B.25) 



where q'2 is given by Eq. (HJi. The secondary miiTor region at (fx 
with radial coordinate between r2 and ^0 is then obstructed: this 
is equivalent to an obscuration on the primary mirror between 
z - Q and some z - Zx < Lx (Fig. |4]i. The correspondence 
between ri and r2 is given by the equation (Spiga et al. 12009b 



r2 ^ 



2{a Q + ax)RQ - 3axrx 



(B.26) 



By comparing Eqs. ( IB.25b and ( IB.26b . we find rx, and the corre- 
sponding value of Zx using Eq. (IB. lb . 



Zx ^ 



(Qf2 + 3^0)^2 ~ (^0 - Rm) 

ax 



(B.27) 



from which, recalling that 02 = 20-0 - ci (Eq. (|5]l), we obtain the 
vignetting coefficient, V3 = 1 -Zx/Lx, 



Viiax) - 1 + 



- a2) 
Lxax 



(B.28) 



where we have used the definition of 2 (Eq. (|25]) ). We have 
thereby found Eq. ([32]). 
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